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^\ ■ The Bethe ansatz equation is solved to obtain analytically the leading finite-size correction of 

the spectra of the asymmetric XXZ chain and the accompanying isotropic 6-vertex model near 
the antiferromagnetic phase boundary at zero vertical field. The energy gaps scale with size N as 
iV" 1//2 and their amplitudes are obtained in terms of level-dependent scaling functions. Exactly on 
the phase boundary, the amplitudes are proportional to a sum of square-root of integers and an 
anomaly term. By summing over all low-lying levels, the partition functions are obtained explicitly. 
• Similar analysis is performed also at the phase boundary of zero horizontal field in which case the 

energy gaps scale as iV -2 . The partition functions for this case are found to be that of a non- 
relativistic free fermion system. From symmetry of the lattice model under n/2 rotation, several 
' identities between the partition functions are found. The iV -1 / 2 scaling at zero vertical field is 

interpreted as a feature arising from viewing the Pokrovsky-Talapov transition with the space and 
time coordinates interchanged. 
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03 ■ I. INTRODUCTION 



Spectral properties of the asymmetric XXZ chain and its associated lattice model, the asymmetric 6-vertex model, 
have been of much interest recently The asymmetric XXZ chain is a non-hermitian generalization of the 

standard spin 1/2 XXZ chain || and is the anisotropic limit of the row-to-row transfer matrix of the asymmetric 
6-vertex model, i.e. the general 6-vertex model in horizontal and vertical fields Both models are solvable by the 

Bethe ansatz for arbitrary interaction parameter A, horizontal field H and vertical field V, and for the lattice model, 
also for arbitrary anisotropy parameter. 

When the interaction parameter A is less that —1, the asymmetric XXZ chain and the asymmetric 6-vertex model 
Q\ • are antiferromagnetically ordered for small values of H and V. (We use the term antiferromagnetic phase to denote the 
\Q \ antiferroelectric phase of the lattice model.) The ordered phase in the H-V plane is bounded by the antiferromagnetic 
• phase boundary beyond which the systems are disordered. The disordered phase is critical in that excitations are 
massless and correlation lengths decay algebraically. For still larger values of fields, complete ferromagnetic ordering 
sets in. The free energy of the 6-vertex model as a function of the two fields can be interpreted as the equilibrium 
crystal shape PJlO| and the antiferromagnetically ordered region of fields corresponds to the flat facet, and the 
antiferromagnetic phase boundary to the facet boundary. On the facet boundary, the curvature of the free energy 
surface is universal Jll|,p|. The nature of phase transitions at the facet boundary can be understood from the view 
point of domain wall excitations and is that of the Pokrovsky-Talapov (PT) transition ]l^ , |l3| |. 
I ■ Within the critical phase and on the phase boundaries, there exist infinite number of massless excitations in the 
thermodynamic limit. When the system size, i.e. the number of spins in the chain N, is finite, the degeneracies are 
lifted producing energy gaps. Finite-size scaling of these energy gaps gives valuable information on the properties of 
the system [Q. Throughout the critical phase of the asymmetric XXZ chain and the asymmetric 6-vertex model, the 
energy gaps scale as iV -1 , and as shown in [0, the finite-size scaling amplitudes can be accounted for by the central 



charge c=l conformal field theory |15|, with suitable modifications to account for incommensurate arrow densities. 
All finite-size scaling information are encoded in the 0(1) part of the partition function which is covariant under the 
modular transformations and, for the case of the asymmetric 6-vertex model, is given by the modified Coulombic 
partition function. 

When A > 1 and V = 0, the critical phase is bounded by the so-called stochastic line where the asymmetric XXZ 
chain Hamiltonian describes time evolution of the single-step model jH| which is one of the simplest realization of the 
Kardar-Parisi-Zhang (KPZ) universality class of non-equilibrium growth models pi] and also of the driven lattice gas 
Jig po| . On the stochastic line, the TV -1 scaling crossovers to the KPZ type N^i 2 scalingwhere the exponent 3/2 
is the dynamic exponent for the 1+1 dimensional KPZ class. This has been shown first in [|3| for a special limit, and 
numerically in j| for other cases. More recently, it was shown in [||, which will be denoted as I in this work, that the 
energy gaps for general cases can be systematically expanded in power series of TV" 1 / 2 and the scaling amplitudes are 
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given by level-dependent, but universal, scaling functions. 

Recently, the finite-size scaling on the antiferromagnetic phase boundary of V = is discussed by Albertini et al. 

0. In their work, it was shown analytically that the low- lying excitations satisfy the energy- momentum relation 
e ~ y/ik in the thermodynamic limit which in turn indicates the A -1 / 2 scaling of energy gaps. That the energy gaps 
scale as TV -1 / 2 was then confirmed by numerical solutions of the Bethe ansatz equations for finite N. 

In this work, we present analytic solutions for the finite-size scaling at the antiferromagnetic phase boundary using 
the method developed in I. It is seen that the A -1 / 2 scaling of energy gaps for V = phase boundary arises from the 
same mathematical origin as in the case of the stochastic line. The scaling amplitudes for general levels are given in 
terms of a scaling variable u ~ (H — H C )N where H c is the position of the phase boundary at V = 0. Their crossover 
behaviors enable one to identify the class of levels which become degenerate at the phase boundary of \H | < H c . At 
H = H Cl the energy gap amplitudes are found to be proportional to a sum of square-root of integers and an anomaly 
term. Furthermore, the O(l) part of the partition function in the scaling limit, which is simply called the partition 
function in this work, is also obtained for both the chain and the lattice model in terms of several infinite products 
involving non-integer powers of nome. Similar analysis is also carried out at the phase boundary of H = 0. Here, the 
energy gaps scale as N~ 2 which is a characteristic of the PT transition. Also the corresponding partition functions 
are evaluated and shown to be that of the non-relativistic free fermion system with dispersion relation e ~ k 2 . Using 
physical requirement of invariance of the partition function of the lattice model under exchange of rows and columns, 

1. e. under tt/2 rotation of the lattice, we find interesting mathematical identities between infinite products. We also 
argue that the A -1 / 2 scaling at the V = phase boundary is a feature arising from viewing the standard PT transition 
of the H = phase boundary after 7r/2 rotation. 

This paper is organized as follows. In Section 0, we set notations and review the Bethe ansatz method for bulk 
properties. Then we discuss classifications of low- lying excitations. In Section [II, the formalism for calculating 
finite-size corrections in energy is developed in line with I and the leading order solution is expressed in terms of 
level-dependent scaling functions with the scaling variable u. In Section IV, the crossover behaviors of the energy 
gaps for H > H c and H < H c are discussed. From this, the class of excitations which remain massless at the facet 
boundary of \H\ < H c is identified. In Section the spectra and the partition functions at H = H c are derived. 
We discuss in Section VI, the spectra and the partition functions at the antiferromagnetic phase boundary of H = 0. 
From t his, w e obtain several identities between infinite products. Finally, we summarize and discuss our results in 
Section VII. Some mathematical details are relegated to Appendices. Appendix ^ proves some properties of the 
scaling function, Appendix ^ shows the details of derivations for the partition function at the V = phase boundary, 
while Appendices |y and [S] derive the spectra and the partition function, respectively, at the H = phase boundary. 



II. BETHE ANSATZ AND ENERGY LEVELS 



We consider the asymmetric XXZ chain of N sites whose Hamiltonian is given by 

H = (2sinh A)- 1 f){£E£* WoJh _ l) _ e™at*- +1 e-^^i} (D 

where af are the Pauli spin operators, cr^ r+1 = cr", and H is the horizontal field. The standard interaction parameter 
A has been parametrized as A = — cosh A with A > since we are interested in the region A < — 1. The front factor 
is included for later conveniences. Below we use the short notation 

a = exp(-2A). (2) 

We also use the notation for the magnetization as 

N 

J>*=A-2Q = 2r (3) 
»=i 

where Q is the number of down spins and use 

g = g/* = i-£ (4) 

so that < q < 1 and r takes the integer (half-integer) values for N even (odd) . We work in sectors of general but 
finite r as N — > oo. Since the spectra of Eq. (0) is symmetric with respect to H and r, we work in the region of H > 
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and r > 0. When a vertical field V is present, it contributes an additional energy — 2rV in Eq. (|l|) but we set V — 
except in Section VI where the spectra of the H — phase boundary are discussed. 
The eigenvalues E of Tt in the sector Q are given by 

Q 

E = J2fE(x,) (5) 

3=1 

where 

/*(*)=— ^-t^-, (6) 
a; — 1 1 — ax 

and {ay} are the solution of the Bethe ansatz equation 

N Q 



1 — axi 



(.ijg-i JJ ±!_^ (i = l,2,...,0). (7) 



3 = 1 J 



(The roots {x-,} are related to the {a-,} notation of by Xj = exp(A — iaj). ) 

The asymmetric XXZ Hamiltonian is the logarithmic derivative of the asymmetric 6-vertex model row-to-row 
transfer matrix at the extreme anisotropic limit jq]. The eigenvalues A of the latter are also expressed in terms of 
{xj}. In this work, we consider the isotropic 6-vertex model, or the F model, in external fields whose six Boltzmann 
weights are given as = exp(±(i? + V)), w 3j 4 = exp(±(if — V)) and ^5,6 = 2 cosh(A/2), respectively. The transfer 
matrix eigenvalues of the F model are given as 

Q 

-lnA= -HN -2rV -^2f F ( Xj ) (8) 

3=1 

where 

/ F (x) = In — — (9) 

1 — CXp( — A)x 

for H > 0. Following discussions can be easily extended to the anisotropic 6-vertex model but we consider only the 
isotropic case for simplicity. We define the phase function or the counting function Zm{x) by 

1 Q 

iZ N (x) = 2H - A + (1 - q) Ins + In + — V / z (g, gj) (10) 

— rvr /V ' ^ 



1 - ax AT • 

3=1 



with 



/zOr, x') = Ins' - In ) ^ = lnx' + £ — f ^ V (11) 

1 — ax x n \ x I 

and the root density function Rn(x) by the derivative of Zn as 

Rn{x) — ixZ' N {x) 

1 Q 

-^/ fl (x,^) (12) 



Q 

x ax 1 

-9 



x — 1 1 — ax A" 

3=1 



where 



W v ^^ Q w(i| (13) 



with X^n^o denoting the sum over all integers except 0. In terms of Zn, the Bethe ansatz equation can be rewritten 



as 
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Z N (x j ) = ^I j (I -I) 



where Ij are half-integers (integers) for Q even (odd). The ground state for each Q is obtained if one chooses Ij as 



lix 2-k ( Q + l 

iv j ~ 77 



*i = TCF H^+i Mi. (15) 



for j=l,2,. . In the thermodynamic limit, {xj} for the ground state form a continuous curve given by the locus 
x = Z^((j)) for — nq < (f> < irq where Z^ 1 (</>) denotes the inverse function of Z N (x). 

The simplifying feature for the case of q = 1/2 and \H\ < H c is that the root contour becomes closed in the x-plane 
as for the case of the stochastic line. Here, H c is the antiferromagnetic transition point to be given in Eq. (^l|). To 
find actual form of Z OD (x) for q = 1/2 and \H\ < H c , one may evaluate the sum in Eq. (|l(i| ) by a contour integration 
over the circle \x\ = exp(A — 6) with — A < b < A using the known expression Roc^x) which is 



M*) I—'''" ( 16 ) 



with J2 n denoting the sum over all integers and 1 < |x| < 1/ar. The result is 

iZ^x) = 2H-X + -lnx + -ln\x \ + J2 n , n d x n + x%) (17) 
2 2 ^„(i + a n) 

where xo = exp(A — b + ni) is the end point of the contour. The conditions Zoo(exp(A — b ± ni)) = ±7r/2 are met if 
H is related to xq or 6 by 

2if = A-lnN-2^— — (18) 

It is convenient to use the real function 5(6) defined by ]7],|| 

m /,s 6 (— l) n sinhn6 

3(6) = - + > - — 

2 ^— ' n cosh nA 

n— 1 

_ ln Cosh(A + 6)/2 6 y, (-l)« a ™sinhn6 

~~ 11 cosh(A - 6)/2 ~ 2 ~ ncoshnA ' ^ ' 

Eq. ( |l8| ) is then the same as 

H = 5(6) (20) 

while i?oo(^o) = 2' (6) with ' denoting the derivative with respect to 6. 

As 6 varies from to A, xq moves from — e A to —1 and increases from to 

H c = 5(A). (21) 

It turns out that the critical value H c is indeed the value of the antiferromagnetic phase boundary at V = 0. Moreover 
the entire antiferromagnetic phase boundary of the lattice model is described by the ii-dependent critical vertical 
field V c given by @ 

V c = S(A - 6) (lattice) (22) 

with — 2A < 6 < 2A. That of the asymmetric XXZ chain with normalization as given in Eq. ([|) is given by (See 
Appendix C) 

V c = 5' (6) (chain). (23) 

Typical phase diagrams in the H-V plane are shown in Fig. 1. 

The bulk ground state energy per site for \H\ < H c is given as eoo = liiriAr^oo E/N = — ep where 
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1 n 

e = -|2V- (24) 

2 ^ 1 + a n v ' 

n—l 

and the bulk free energy of the F model is /oo = limjv^ co (— lnA)/JV = —A/2 — Y^=i ™~ 1 a Tl/ ' 2 tanhnA. 

In the antiferromagnetic phase, finite mass gaps appear in the energy spectra. But these mass gaps are offset 
by fields and at the antiferromagnetic phase boundary, infinite number of other levels become degenerate with the 
ground state in the thermodynamic limit. These degeneracies are lifted for finite but large N. Our aim here is to find 
these finite-size corrections for arbitrary energy levels. We denote these energy gaps as AE and AF, respectively: 
AE = E - e^N and AF = — In A — f^N . 

The ground state characterized by Eq. ( |l5| ) corresponds to a Fermi sea as depicted in Fig. 2(a) for r=0. Other 
energy levels are obtained if other sets of are chosen. An important class of levels is the m-shifted levels where 
{Ij} are shifted by an integer m from that of the ground state, i.e. 2irIj/N — 4>j+ m - These states are denoted as 
(r,m,0). Fig. 2(b) shows the 1-shifted level (0,1,0). Further excitations are obtained by creating from the m-shifted 
states an equal numbers of particles and holes near the two ends of the Fermi sea. Excited levels are then characterized 
by their particle and hole positions. The particle (hole) positions near the right-hand side of the Fermi sea can be 
specified by a set of integers {pk}, {{hk}) which are related to the index j in <fij as j = Q + m + pk (Q + m + 1 — hk), 
for fc=l,2,. . n + being the number of particle-hole pairs, and 1 < p\ < p2 < ■ ■ ■ < p n+ (1 < h\ < h% < ■ ■ ■ < h n+ ). 
Particles (holes) near the left-hand side may be labeled similarly by {pk{ ( {hk}) with j = 1 + m — pk (m + hk), 
for fc=l,2,. . .,n_ and 1 < p\ < p2 < ■ ■ ■ < p n _ (1 < hi < hi < • • • < h n _) . Such excitation configurations are 
denoted collectively by V . General levels are then labeled as (r,m, V). Fig. 2(c) shows a level obtained by creating 
n + =2 particle- hole pairs at the right-end and n_=l pair at the left-end from the 1-shifted level of Fig. 2(b). In Q, it 
was shown for the critical region that levels obtained in this way account for all the low-lying excitations which are 
expected to appear from the central charge c = 1 conformal field theory. Here, we also assume that all the low-lying 
excitations are generated in this manner. 

When the particle-hole configurations differ from the ground state only at finite distances away from the two ends 
of the Fermi sea, finite-size corrections in energies arc determined by the analytic property of ^00(2;) near x = xq. 
However, the root density vanishes at x — x° c = e 7 ™ which is the end point of the root contour for H = H c . Therefore 
Z^((f>) exhibits a square root singularity at <j> — ±n/2. The same phenomenon was the origin of the unusual 7V~ 3 / 2 
scaling of the energy gaps along the stochastic line and one can expect similar mechanism produces N~ x / 2 scaling 
in the present case. On the other hand, for \H\ < H c , 0(N~ 1 ) variations in </> near <f> — ±tt/2 cause, via Z^ 1 (4>), 
0(N~ 1 ) variations in x-plane and hence the finite-size corrections take the form of power series in A^ 1 . What is 
special at the antiferromagnetic phase boundary is that the real parts of the energy gaps scale as Af -2 which is the 
characteristic of the PT transition. These features are born out explicitly below. 



III. ENERGY GAPS FOR H NEAR H c 



In this section, we derive the leading order finite-size correction of the energy gaps AE and AF for general levels 



with (H — H C )N as a scaling variable. Readers not interested in details may skip to Eqs. ([|2|) and (44) which are the 
main results of this section. 

To find the energy gaps for H near H c , we employ the method of I. Thus, for finite but large N and for levels 
whose energies are close to that of the ground state, we assume, following I, that there exists in the complex x-plane 
a point x c near x° c = e 7 ™ such that Z' N (x c ) = 0, Zjy(x) itself having a branch cut passing through x c , and that near 
x c , Zn(x) has the expansion of the form 



iZ N (x) = ±iirq + y±^ + ^{x - x c f + 0((x - x c ) 3 ) (25) 
A af 



7T 1 

A 

where the upper (lower) sign refers to the assumed value of iZjq(x) at x = x c (x = x c e~ 27rz ). This amounts to 
assuming 

Z N \± 7 r q + ^ = xt± l a 1 ^^^+0(N- 1 ) (26) 

for £ — O(l) with x+ = x c and x~ = x c e~ 2m . 

The constants x c , a 1; y±, etc. are all level dependent and should be determined self-consistently. But, first we 
assume they are known and evaluate the finite sums of the form 



5 



Q Q 2 

S[f] = E f( x i) = E f{ Z N\^h)) (27) 

2=1 j=l 

to 0(1/">/N) for /(a;') given in Eqs. (§),(§, (0) and (||). Following the same steps as in I, but generalizing to levels 
with r7i 7^ 0, we find 

S[f] = /to^tof + ( m + - (™ + ^-)/(^e" 27ri ) 

+ai/'(x c )r 1 ( 2 / + ,y_)^+0(iV- 1 ) (28) 

where the integral is over a closed contour in the annulus 1 < \x\ < 1/a and Y\(y+, y_) is defined below. To define Y\ 
, we first introduce two functions J+{y) and J-{y) of complex variable defined for Qy > and 9?/ < 0, respectively, 
as 

with the branch cut of the square roots at the negative real axis. As shown in Appendix [A], J± (y) satisfy the recursion 
relations 

J±(y±2i) = J±(y)TH/y- (30) 

Thus, we may extend the definition of J+(y) (J-(y)) into the half plane 3y < (Sy > 0) using the relation Eq. ( |30| ) 
recursively. Having defined J±(y), we then give the expression of Y\ for general level (r, m,V): 



Y\ (y + , y_ ) = J + + i — 27m) + J_ (y_ — i — 2im) 

">■+ 
-i 1 
fe=i 



;E + i(i - 2m - 2 Pfe ) - Vy+ - K 1 + 2m - 2h K 



-i 

k=l 



E y/v- -i(l + 2m-2p fc )- + i(l - 2m - 2fe fc ) . (31) 



The first two terms in Eq. ([n]) are contributions for the m-shifted levels and the sums in the second (third) line 
accounts for particle-hole pairs created at the right- (left-) end of the Fermi sea. When m = and y + =y~, Yi(y,y) 
reduces to minus of Y\(y) defined in I. 

Applying Eq. ( pSj ) to Eq. (12), we then obtain the solution for Rn as 

«»W = *.W + £ + i£ ^ ( J)" + i^.^/I) + o (W -). (32, 

Having obtained i?jv(x) to 0(A~ 3 / 2 ), we can evaluate to 0(N 1 I 2 ) for any f(x). In particular, applying the 
general sum formula to f(x') = fz(x, %'), Zn{x) can be evaluated to Q(N~ 3 / 2 ). Evaluating the latter at x — x c and 
x = x c exp(— 27T7), we obtain two self-consistency equations (See Eq. (|25|).) 

± t7r« + ^y± = ±t(l - g)7r + 2H - A + 2 E J 



- + a n n 

+(1 H lna; r — 7n) + (2im 1 — 

K N 2Ni ' v y v 2 ^iV 

+27r 1 / 2 e a 1 ^ 1 Y 1 ( 2/+ , y_)A- 3 / 2 + (9(A- 2 ) (33) 



where eo is given in Eq. (|2J). From Eq. (|33j), one first obtains y+ — y_ = 4ir. Next, evaluating Rn(x) and R'^{x) at 
x c , we obtain two more self-consistency equations as 

= i?co(£ c ) + ^ + O(A~ 2 ), (34) 
^=i4(*c) + 0(A- 3 / 2 ) (35) 



G 



where y+ — y_ = 4ir has been used in Eq. (|34|). Solving Eqs. (|34|) and ( J35| ) perturbatively, we have 



e? A 

V2 , 



O(A^), (36) 



ai = — + 0(iV- 1 ) (37) 
ei 



where e\ is a positive real constant defined as 

ei ee (-iC(-l)) 1/a = (-2" (A)) 1//2 = f i + 2^ n(-l)" I ^ r J . (38) 

Using Eqs. fl36| ) and ( p7| ) in Eq. p3]), we find the leading order solution of y± as 

y± = (±2r + m)i+ u + (9(iV- 1/2 ) (39) 

where u is the scaling variable 

u=(H - H c )N/n (40) 

which is assumed fixed as A — > oo. Having determined y±, x c and ai to necessary orders, we are now in a position 
to evaluate Ai? and AF. Applying Eq. @ together with Eq. @ to Eq. (S), we find 



p , y+ -y- 

E = -e A -H 

min(l, a - ™) 



+5=^^' +o(jv "»- <4i) 

Inserting the solutions Eqs. (|3^), ( |37j ) and (|39|), to Eq. (f4l|), the O(l) terms in Eq. ( |4l| ) cancel out and we finally 
have, for V = 0, 

AS ee £ - 6oo iV = ex y[ r ' m ' r) ^m+ O(N^) (42) 

where eoo = — eo, ei is given in Eq. © and y [ {> m < v) is the value of Yi(y+, yJ) at y± = u + (±2r + m)i: 
y(r,m,v) = j + {y t + i{2r- m + l)) + J-{u-i{2r + m+l)) 

{\/ u + *( 2r -m+1- 2p k ) - y/u + i(2r -m-l + 2h k 

fc=i 



k=l 



^2[V U - + m + l- 2p k ) -yju- i(2r + m - 1 + 2h k ) ) . (43) 



Here, J± (y) and it are given in Eq. @ and @, respectively. j;^'™'^ regarded as a function of the scaling variable 
u is then the scaling function. Similar calculation for the F model leads, for V = 0, 

AF ee — In A — f^N = -mm + e _ly^ m > v \\^L + OW" 1 ) (44) 

ei V A" 



where the constant is given as 



1 f— 11" 

n— 1 



Note that AS has an imaginary term due to A being negative for m odd. Apart from this, all energy gaps scale as 
A -1 / 2 and all the level dependences are encoded in the scaling function y( r ' m:V \ Higher order corrections which 
come as a power series in A -1 / 2 can be calculated perturbatively using the method similar to that of I. 
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IV. CROSSOVER BEHAVIORS FOR H =fi He 



Eqs. (^2|) and ([Ij) are obtained with the scaling variable u fixed. By considering the limits u — > ±00, we can derive 
how the spectra crossover to those of the critical and massive phases. First consider the case H > H c where u — > 00 
as N — » 00. From Eq. (|29|), one can show that, as u — > 00, 

J±(« ± to) = -iu 3 / 2 1)^/2 + _ 2 + _L )u -V 2 + 0(u -3/ 2) . (46) 

Using this in Eq. (El), 3;( r,m '' p ) becomes, 

y[ r ' m ' V) = —u 3 ' 2 + imu 1 ' 2 + + r 2 - JL + + AT)^ 1 / 2 + 0( U - 3 / 2 ) (47) 

with Af = J2k=i(Pk + hk — 1) and TV* = X)fc=i(Pfc + ^fc — !)■ The first term in Eq. ( ]47| ) contributes a bulk energy 
term proportional to (F — H C ) 3 ^ 2 N to AF and AF, the second an O(l) imaginary term and the third a real one 
proportional to (F — F c ) -1 / 2 TV -1 . These are exactly those which appear in the critical phase as shown in Q and 
show the complete operator content of the Gaussian model. In the Coulomb gas picture, m and r are the spin wave 
and vortex quantum numbers, respectively, characterizing primary operators of the c = 1 conformal field theory 
whose Gaussian coupling constant g is 2, and JV and J7 account for the conformal towers generated from the primary 
operators. Higher order corrections in Eqs. ([l2|) and (EJJ) not discussed here would give terms which are higher orders 
in (F-F c ). 

Next consider the case H < H c . Now u — ooasiV^oo and one has to be careful about the branch cut of the 
square root. We find 

J±(u ± to) = ±^M 3/2 + G±M 1/2 + 0(M~ 1/2 ) (48) 

with G± = \a — l|/2 for a > or for a negative odd integers but, for a < and even integers, G+ = (|a| — l)/2 while 
G_ = (|a| + 3)/2. Using this, we find the leading order of y^ m,v ^ for levels (r,m,0) as 3^ r ' m '°) = max(2r, ImDIul 1 / 2 
if (2r + m) is even and y^ r ' m ^ = max(2r, \m— 1 1 ) if (2r + m) is odd. Thus, when m is in the range — (2r + 1) < 
m < (2r + 1), yf"™' ) = 2r|u| 1 / 2 + 0(|u| -1 / 2 ). As to the contributions from particle-hole excitations in Eq. pSj), 
we note that, for — (2r + 1) < m < (2r + 1), particles at positions satisfying the conditions 2r — m + 1 — 2pk > or 
2r + m+ 1 — 2j5fe > contribute — lul 1 / 2 but those violating the conditions and all holes contribute | | x / 2 . Thus when 
all particle positions satisfy the conditions, or equivalently, when the last positions satisfy the condition 

- (2r + 1) < m < (2r + l), 2p n+ < 2r - m + I, 2p n _ < 2r + m + I, (49) 

then all particle contributions are cancelled by holes. Thus, for such levels, y^'™ 1 '^ = 2r|u| 1 / 2 + 0(|u| -1 / 2 ) and AF 
becomes 

AF = 2v / 2ei (H c - H) 1/2 r (50) 

as H — > F^ c in the thermodynamic limit. This is the mass gap of the antiferromagnetic phase of the asymmetric XXZ 
chain. The mass gap for general H is derive d in jlj] for levels (r, 0, 0) and is in agreement with Eq. (J50|) . It is also 



derived for general H implicitly in Eq. (C12) and is in accord with Eq. (|23|). Corresponding mass gap for the lattice 
model is 

AF = -i7rm + 2V2e 2 e7 1 (F c - H) 1/2 r. (51) 

For levels which do not satisfy Eq. (|4g|), the mass gap becomes higher by 0(1). Eqs. ( |5C)| ) and (|5|) are valid for r > 0; 
they are valid for general r if r is replaced by \r\. When the vertical field V is applied, both the chain and lattice 
model energy acquire a trivial term — 2rU. When 2rU is equal to the mass gap given in Eq. (|5^) and ( pl|) for the 
chain and lattice model, respectively, all levels (r,m,V) satisfying Eq. (Ej|) and r > become degenerate with the 
ground state energy in the thermodynamic limit. Finite-size corrections for these states will be discussed in Section 



VI 



A special mention is needed for r = which is the half-filled sector for N even. Here, only two levels (0, 0, 0) and 
(0, 1, 0) satisfy Eq. ((49|). The latter which has momemtum ir is the exponentially degenerate second largest eigenvalue 
of the transfer matrix as discussed in p2| . Such exponential degeneracy is reflected in the large u behavour of y^ r,m '' p ^ . 
When expanded as a power series in both j/^ ' 0,0 ' an d yt ' 1 - ) have vanishing coefficients to all orders in the 

series. 
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V. SPECTRA AND PARTITION FUNCTIONS AT H c 



Now we consider the case H = H c (The points A in Fig. 1). Since u = in this case, one needs to evaluate J±(ia) 
for integer values of a. Successive application of the recursion relation Eq. (p0|) yields, 



i E^=o 3)/2 V 1 + 2 i - v^ci if a > and odd 



T , , . /^E£o 1)72 71^27 - V=fci if a < and odd 

J+M = ^ ' ( _2)/2 /jtt , /— r . f ^ n , ( 52 ) 

v~ «z^j=o v2j + v-iC2 it a > and even 
%/! X) j=o 2 v^J + \/~* c 2 if a < and even 

and J—(— ia) = J+(ia)* where the positive numbers ci and C2, called anomaly, are given as 

Cl = -J +il )/V^=---J o dt (53) 

c 2 = J + (0)/V = i = C(3/2)/(2V27r) = 0.293 995 52 .. . (54) 

with £(z) the Riemann zeta function. We will see later in an indirect way that c\ = (v2 — 1)C(3/2)/(4tt) = 
0.086109 29... (Eq. @). 

Using Eq. j52] ) in Eq. (^3|) with it = 0, one can read off y^ r,m ' V ^ for arbitrary levels. It is a sum of square-root of 
integers and the anomaly term. We give a few examples: 
(i) Ground state in each sector r for N even: 

y[ r '°' 0) = V2 fc^TTTj - cA . (55) 



(ii) m-shifted levels in sector r = 0: 



,,(o,m,o) _ f 2-\/i(l + v3 + • • • + \Jm— 1) — v^ci for 771 > even , . 

^ ~~ \ 2y/i(y/2 + + + V2c 2 for to > odd. 1 j 



For m < 0, j;' '" 1,0 ) i s the complex conjugate of yj°'l m l>°). 
(iii) m-shifted levels with — (2r + 1) < to < (2r + 1): 



r-m/2-1 r+TO/2-1 
3; Cnm,o) = _ V 2 Ci+V=I £ VT+27+^ E (57) 



for (2r + to) even and 



r-(m+l)/2 r+(m-l)/2 



y (r,rn,0 )=V 2 c2 + y=I y y27 +V " y 72j (58) 



for (2r + to) odd. 

(iv) One particle-hole pair excitation from (0,0,0): 

,(0,0,P) 



y (o,o,v) = _^ ci + Jl^j 2pi _ I + V^v/2/n - 1 (59) 

with the positive integers p\ (h\) being the position of the particle (hole). 

We have compared above predictions with existing numerical data of AE obtained by solving the Bethe ansatz 
equation for finite N Extrapolations of available data (N < 80 for levels (0,0,0) and (1,0,0), N < 68 for (0,2,0), 

N < 56 for (2,0,0)) all show excellent convergence to the expected exact values. 

Having obtained the leading finite-size corrections for arbitrary levels, one can go further to evaluate the partition 
functions in the scaling limit. We consider the partition function of the asymmetric XXZ chain at temperature T 
(in units of the Boltzmann constant) and the lattice model on N x M lattice with periodic boundary conditions 
defined by Ze = Yl ex P( — AE/T) and Zp = J^exp(— MAF), respectively, where the sums are over all states and 







the bulk part exp(— Neoo/T) for the chain and exp(— MNfoo) for the lattice model has been taken out. To obtain a 
non-trivial limit, we take the thermodynamic limit with \/~NT or MN~ X / 2 fixed. In this scaling limit, the 0(A~ 1//2 ) 
part of the energy gaps produces finite contributions while levels which are not degenerate with the ground state in 
the thermodynamic limit need not be considered. 

Considering Zp first, and we can write from Eq. (E4J), 



Z F = ^exp(-MAF) 

£ (-l)^exp(-^yfy^ 



(60) 



(r,m,P) 

The sum is conveniently expressed in terms of a complex parameter q which we call nomc: 

q = cxp(— nVir) (61) 

where 

t = 2e 2 2 M 2 /(-Ke\N). (62) 

Derivation of the partition function is relegated to Appendix |b[ Our final result for Zp is 

/ Zi(q) ± Z 2 (q) + Z 3 (q) for M even 
\ ±Z 1 (q)+Z 2 (q)T2 3 (q) for M odd 

where the upper (lower) sign is for N even (odd) and Zj(q) are defined as 



(63) 



OO 

Z x (q) = -(qq)" ci J] (l + q 71 ^) (l + q 7 ^) (64) 

J=-0O 

Z 2 (q) = l(qq)- c i (l - q 7 ^) (l - q 7 ^) (65) 

j=— OO 
OO 

2 3 (q) = ^(qq) C2 II f 1 + ^) i 1 + ^) ( 66 ) 



with c\ and c 2 given in Eqs. ( p3| ) and (|54|), respectively, and q denoting the complex conjugate of q. 
For the asymmetric XXZ chain, we have 

Ze= E expC-l^/fM^)- (67) 

In this case, we can simply use the result of Zp for M even by redefining the nome as 

, ei /2tti 

q =cxp (~rV^' ^ 

Using this, we then have 

Z B = Z 1 (q')±Z 2 (q') + Z 3 (q') (69) 
where the -I- (— ) sign is for N even (odd). 



VI. SPECTRA AND PARTITION FUNCTIONS AT H = AND V = V c 



The partition function of the isotropic lattice model on N x M lattice should be invariant under rotation of the 
lattice by 90°. After the rotation, roles of H and V are also interchanged. Thus we expect an invariance of the partition 
functions under N <-» M and H «-> V. When the system is in the critical phase, this symmetry is manifest by the 
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modular covariance of the partition function. Modular invariance plays an important role in classifying the conformal 
invariant theories and in revealing their mathematical structures. Present system is in a way a non-relativistic limit 
of the CFT and it is of interest to find the analog of the modular transformation. For this purpose, we have also 
calculated the energy spectra of the two models at \H\ < H c for those which are low- lying near V = V C where V c is 
the ii-dependent critical vertical field. The method employed is similar to the case of H = H c and V = except that 
the ansatz for Z^ 1 (</>) takes a different form and that only those levels satisfying the condition Eq. ([i"9| ) and r > 
come into the working. Details of the derivations are given in Appendix ^|. To compare the result with that of the 
V = phase boundary (the points A in Fig. 1), we consider only the H = phase boundary (the points B in Fig. 1). 
Our result for the F model at H — is 



AF = -i:im - 2r(V - H c ) + ' m ' V) (n/N) 2 + 0(N~ 3 ) 

2e 2 

where H c , e\ and e 2 are given in Eqs. (pl|), (38) and (fi"5|), respectively, and y^'™ 1 '^ is given by 



(70) 



y{r,m,V) 



2rm 



E 

E 

fc=i 



((2r 
((2r- 



while that for the chain at H = is 

AE = 2re 2 



1 - 2 Pk 



1 - 2p k 



64 , 7 (r.m,-p) 

2?/ 4 



(2r 
(2r 



2h k ) 2 ) 



l + 2hu 



(tt/N) 2 + 0(N~ 3 ) 



(71) 



(72) 



where = — S"'(0) = J2^=i(~ l) n+1 n 2 / coshnA. Note that V c at H = is the same as H c for the lattice due to 
isotropy, but it is e 2 = S'(0) for the chain. For H not equal to zero, imaginary 0(N~ 1 ) terms appear in AE and AF. 
However, the real parts are 0(N~ 2 ) for all \H\ < H c . 

From the result of the energy gaps, we are able to calculate the partition functions as in the previous section. We 
consider the lattice case first since it is more general. To compare with Eq. (|63|), we interchange M and N in Eq. j70| ) 
and offset the mass gap by setting V = H c . Then the partition function of a M x N lattice at the phase boundary is 



Z F = 



E 

(r,m,'P) 



(-1) 



mN 



CX P ( 3^4 ) 

T 



(73) 



with t — 2e 2 M 2 /(ire 2 N) as was given previously in Eq. (|(32|). This quantity should be the same as that defined in 
Eq. ( |60| ) by symmetry. To express the partition function in a compact way, we define the new nome as 



exp(— 7t/t). 



(74) 



This is analogous to the nome corresponding to the conjugate modulus in elliptic functions |6). After performing the 
sums over the levels, we find Zp as 



Zi(p) ± Z 2 (p) + Z 3 {p) for M even 
±Z 1 (p) + Z 2 {p) T Z 3 (p) for M odd 



(75) 



where the upper (lower) sign is for N even (odd) and Zi are defined as 



Zi(p) 
Z 2 (p) 

Mp) 



1 oo 

i nc 

j=-oo 
1 oo 

5 no 

j=-oo 



,(l+2j) 2 



,(l+2j) 2 



(76) 
(77) 
(78) 
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Derivation of Eq. (f75|) is sketched in Appendix The corresponding partition function Ze for the asymmetric 
XXZ chain with M sites at temperature T can be obtained from the N even result of Eq. ( |75|) with p replaced by 
p' = exp(-e i TT 2 /(2elM 2 T)) ; 

Z E = Z 1 (p')±Z 3 (p')+Z 2 (p') (79) 

with the + (— ) s ign for M even (odd). 

Equating Eq. (j63|) and d73) for the 4 cases of parity of M and AT, we find 



Zi(q) = Zi(p) (80) 

for all i=l, 2, 3, and for any r where q and p are related to r by Eq. ( |6l| ) and ([74]), respectively. This is an analog 
of the conjugate modulus transformation of the elliptic functions. These identities are obtained by calculating the 
physically same quantity in two independent ways and are confirmed numerically. However we are not able to prove 
them directly. A byproduct of the identities is the analytic expression of c\ mentioned below Eq. (154). Taking the 



r — > oo limit in ln^i(exp(— 7r/r)), and evaluating the leading order contributions using Eq. (A3), one easily obtains 

ln^expt-Tr/r)) — ±(V2- l)C(3/2) v^A (81) 



Comparing this with Zi(exp(— 7rvzr)), one obtains the alternate result for ci as 

c 1 = (v^-1)C(3/2)/(4tt). (82) 
Similar steps using Z3(exp(— tt/t)) confirm the leading order behavior of Zs(q) for small q. 



VII. DISCUSSIONS 



In this paper, we have presented Bethe ansatz solutions for the finite-size corrections of the energy gaps at the 
antiferromagnetic phase transition for both the asymmetric XXZ chain and the isotropic 6-vertex model or the F 
model. Furthermore, all the low-lying levels are summed to obtain explicit expressions for partition functions. As a 
byproduct, interesting identities between several infinite products are found. Main results of this paper are Eqs. ( [42] ) 
and (|44|) for energy gaps at V = phase boundary, Eqs. (|t^) and (|70| ) for the same at H = phase boundary, Eqs. 
( |63| ) and ( |7q ) for the partition functions and the identities Eq. (|30|). 

The finite-size corrections of the energy gaps are calculated using the method developed in I ( ||). Here, the phase 
function appearing in the Bethe ansatz solution is assumed to take certain form and is determined self-consistently. 
At H — H c and V = 0, the finite-size corrections take the form of power series in N~ x / 2 because of the fact that 
the density of roots for the ground state vanishes at the end points of the root contour. The same is true on the 
stochastic line treated in I but, in the latter case, the leading 0(N~ 1 ^ 2 ) term is absent resulting in the N~ 3 / 2 scaling 
of energy gaps. For \H\ < H c and V = V C , the finite-size corrections come out as a power series in AT -1 . The leading 
0(N~ 1 ) term is imaginary and the real contribution appears from 0(N~ 2 ). This A^ 2 scaling is the characteristic 
of the PT transition and is more generic. The points H = ±H C and V = of the phase boundary are special in 
that the direction of the 'time' of the Hamiltonian Eq. ([[]) and the row-to-row transfer matrix of the lattice model is 
orthogonal to the direction of the field and tangential to the facet boundary in the H-V plane. 

The partition function at H = and V = V c given in Eq. ( f75| ) is actually that of an one-dimensional non-relativistic 
free fermions with dispersion relation e = Ak 2 , A being a constant. To make exact correspondence with Eq. (f75|), we 
simply require the momenta k be even or odd integer multiple of tt/L, L being the length of the system, depending 
on the parity of the total number of fermions. Thus the effective theory on the H — antiferromagnetic phase 
boundary is that of non-relativistic free fermions with the inverse temperature or the number of rows playing the role 
of imaginary time. This is also true at q = or q = 1 (the dashed lines in Fig. 1) where the transition can be viewed 
as the PT-type commensurate-incommensurate transition ]l3| , p4"| | . For such effective theory, the single particle wave 
function is exp(ikx — et) where t is the imaginary time. Now we consider the effect of rotation in x-t plane. Such 
transformation is not covered in the Schrodinger group |2J| but makes sense in the context of lattice models. When 
the coordinates are rotated by an angle 9, in order for the wave function to remain invariant, the dispersion relation of 
the rotated system must change to e' = ik' tan# + A(cos0)~ 3 k' 2 + 0(fc' 3 ) when 6 ^ ±n/2 but to e' = y / ±ik'/A when 
6 = ±7r/2. Since the lattice model at H = H c and V — is obtained by rotating that at H — and V = V c by n/2, 
one naturally has excitations with yik' dispersion at the V — phase boundary. From this, we can interprete the 
A -1 / 2 scaling at H = H c and V = as the feature arising from viewing the ordinary PT transition with space-time 
coordinate interchanged. However, the anomalies c\ and C2 in the rotated system are not explained by such simple 
effective theory. 
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APPENDIX A: RESURSION RELATION OF J± 

To show the recursion relation Eq. (EM), we introduce an auxiliary function defined for > as 



J(y) 



i 

2i 



y/y + it - y/y - it _ 1 1/2 1 3/2 



(Al) 



J(y) is related to J±(y) by J+(y) = \J~^iJ{—iy) for 3y > and J_(y) = \/iJ(iy) for 3y < 0. Eq. ( p0| ) follows from 
the recursion relation for J(y) which is 



J(y + 2) = J(y) + 



for 5Ry > 0. To show Eq. (A2), we employ the general sum formula 



£/(?) = / 

3=0 J ° 



f(n,t)-f(0,t) 



dt 



(A2) 



(A3) 



where /(a;, <) = (f(x+it) — f(x—it))/2i. This formula is valid if f(x) is analytic and satisfies limt^±oo e 2 ^f(x+it) = 
in the strip —S < $tx < n + 5 for some 8 > 0. Applying Eq. (A3) to a trivial sum X^=o Vv + 2 j , one then obtains 



APPENDIX B: PARTITION FUNCTION AT H — He AND V = 



In this appendix, we evaluate the partition function of the lattice model at the V = phase boundary. That of the 
asymmetric XXZ chain follows simply redefining the nome and using the result of the lattice model for M even. For 
simplicity, we consider the case of M even and mention the other case at the end. Since M is even, from Eq. (|60|), we 
are to evaluate 

z F = cM-^yl r - m ' v) ), (bi) 

where r is given in Eq. (|62]). 

First consider the case (2r + m) even. Further, if m is in the range — (2r + 1) < m < (2r + 1), y l ^ r ' m ^ i s given in 
Eq. @ and we have 

r-m/2-l r+m/2-1 

exp(-W^ m ' 0) ) = (qqri J] ^ II ^ (B2) 

3=0 j=0 

where q = exp(— irs/ir) and q = exp(— it\f~ vr). From Eq. (^), one notes that exp(— TTy/ry < f' m '' P ' > ) is a product 
of exp(— Try/ry[ r ' m : °^) and factors coming from particle and hole excitations. A particle at position pk with 2pk < 
(2r — m + 1) contributes a factor q-v' 1 + 2 J where j = r — m/2 — pk = 0, 1, . . . , (r — m/2 — 1), while that with 
2pk > (2r — m + 1) a factor qv' 1 + 2 J where j = Pk — r + ti/2 — 1 = 0, 1, ... , oo. A hole at /i/c gives qv T + 2 ? where 
j = hk+r — m/2 — 1 = r — m/2, r — m/2 + 1, . . . , oo. Summing over all particle- hole configurations at the right-end of 
the Fermi sea corresponds to evaluating a free fermion grand-canonical partition function under the constraint that 
the total number of particles should be the same as that of holes. If we denote the fugacity of each particle (hole) 
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by z (z 1 ), the desired sum is obtained by projecting out the z° term of the free fermion grand-canonical partition 
function which is 

r— m/2— 1 oo oo 

E 1 {z)= J] (l + ^- VTT71 )U( 1 + ^^ II (l + a-V 1 ^)- (B3) 

3=0 3=0 j=i — m/2 

S„(z) of this Appendix and Appendix [D] are not to be confused with S(6) defined in Eq. ([l9|). We introduce a function 
A{z) by 

OO OO 

A(z) = H(l + z q ^)l[(l + z- 1 c i ^), (B4) 

3=0 3=0 

and denote its Fourier coefficients by A n ; A{z) = J2 n A n z n . Then multipying and dividing by a factor Ilj=o 1 ^ 2 X (l + 
z -iqVT+2j) on t ne right-hand side of Eq. (p33|), we can put Si in the form 

r— m/2 — 1 

Hi(z) = z r - m/2 A(z) JJ q-^+^J. (B5) 

3=0 

Thus the coefficient of z° of Hi(z) is A_ r+m / 2 n^Z^ 1 ^ 2- 1 q _v/1+2j '. Note that the product over j in Eq. (B5) is 
cancelled by that in Eq. (B2). Summation over particle- hole configurations at the left-end of the Fermi sea proceeds 
in exactly the same way with m replaced by —m and q by q. Thus, the partition function after summation over all 
excitations from each m-shifted state is given as 

J2^v(-^y[ r ' m ' v) ) = (mr ci A_ r+m/2 A_ r _ m/2 (B6) 

v 

for — (2r + 1) < m < (2r + 1) where A n stands for the complex conjugate of A n . Similar calculations for each case of 
m > (2r + 1) and m < — (2r + 1) show that Eq. ( p6| ) holds for all m, provided (2r + m) is even. We then sum the 
right-hand side of Eq. (|B6| ) over all even (odd) m's for 2r even (odd). However, we note that ■ / ^-r+m/2^-r—m/2 
is simply the coefficient of z~ 2r in |A(z)| 2 . To obtain the partition function in the (2r + m) even sector, we finally 
sum over all integer (half integer) values of r if N is even (odd). Even though above derivations assumed r > 0, we 
may sum over negative r's by the symmetry, |A(z)| 2 = |A(z _1 )| 2 . Since the sum over r is equivalent to summing over 
all coefficients of even (odd) powers of z in |A(z)| 2 , it can be written as (^4(1) 2 + A(-l) 2 )/2 ((A(l) 2 - A(-l) 2 )/2) for 
N even (odd). Noting that (qq) _Cl A(l) 2 /2 = Zi(q) and (qq) _c M(— 1) 2 /2 = 2 2 (q), we have the partition function in 
the (2r -I- m) even sector as i?i(q) ± Z 2 (q) with + (— ) sign for N even (odd). 

Next consider the case of (2r + m) odd. If m is in the range — (2r + 1) < m < (2r + 1), y^'" 1 ' ^ is given in Eq. ( |5S| ) 
and we have 

r-(m+l)/2 r+(m-l)/2 

exp(-^M' %m ' 0) ) = (qq) C2 J[ ^ J[ ^ ( B7 ) 

3=0 3=0 

A particle at position pk with 2pk < (2r — m + 1) contributes a factor q _v/2 ^ where j = r — (m — l)/2 — pk = 
0, 1, . . . , r — (m + l)/2, while that with 2p^ > (2r — m + 1) a factor q^> where j = pk — r + (m — l)/2 = 1, 2, . . . , oo. 
A hole at hk gives where j = hk +r — (m + 1)/2 = r — (m — l)/2, r — (m — l)/2 + 1, . . . oo. Thus the free fermion 
grand-canonical partition function for excitations at the right-hand side is 

r— (m+l)/2 oo oo 

S2(*)= n (i + ^q" v ^)II( 1 + zc ' v ^) II (i + ^'q^ 7 )- (B8) 

3=0 3=1 j'=r-(m-l)/2 

In line with the (2r + m) even sector, we introduce a function B(z) by 

oo oo 

B(z) = H(l + z^)l[(l + z- 1 ^), (B9) 

3 = 1 3=0 
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and denote its Fourier coficients by B n ; B(z) — J2 n B n z n - Then S2 becomes 



E 2 (z) = z r ^ m -^ 2 B(z) 



r-(m+l)/2 

n < 

3=0 



(BIO) 



and one sees that contribution from particle-hole excitations at the right-end is proportional to -B- r +(m-i)/2- Adding 
that of the left-end, we have the intermediate partition function as 



B_ 



r-(ra+l)/2 



(Bll) 



for — (2r + 1) < to < (2r + 1) and (2r + m) odd. Similar calculations for each case of to > (2r + 1) and m < — (2r + 1) 
agai n sho w that Eq. (Bll) holds for all to, provided (2r + to) is odd. Proceeding as above, we sum over to and r in 
Eq. ( Bll ) using the symmetry z|£?(z)| 2 = z \B(z )\ and finally find the partition function in (2r + m) odd sector 



as (qq) C2 (5(l) 2 + 5(-l) 2 )/2 for 2r odd and (qq) Ca (B{1) 2 - B(-l) 2 )/2 for 2r even. But B{-1) = so that it becomes 



(qq) C2 -B(l) /2 = Z 3 (q) for both parity of N. Adding this to that of (2r + to) even sector, we have Eq. 
even. 



for M 



When M is odd, we need to insert a factor (— l) m to Eqs. (B6) and (Bll). This does not complicate the sums since 
the sums over to in each sector a re in steps of 2 and the sign factor is equivalent to (— l) 2r for the case of Eq. (jB6| ) 
and (— l) 2r+1 for the case of Eq. (Bll). These extra sign factors then give the full partition function as given in Eq. 
(||) for M odd. 



APPENDIX C: SPECTRA AT H = AND V = Vc 

In this section, we de rive AE andAi* 1 at the phase boundary of H = 0. Although our main interest is for H = 0, 
discussions up to Eqs. ( C12| ) and (C13]hold for general < H < H c . 

When \H\ < H c , given by Eq71l7h has non- vanishing first derivative at xq = exp(-7rz + A — 6), the end point of 
the root distribution in bulk limit, which is related to H by Eq. (pfj|). Thus our starting ansatz for Z^ N 1 (x) is, instead 
of Eq. (|2F 



Z- 1 (±nq + n£/N) = 4 + a 2 n(y ± - i£)/N + a A n 2 (y ± - itf/N 2 + ■■■ 



(CI) 



x e 



-2iri 



Here, y±, 0,2 and 04 are to be determined self-consistently. We then follow the 



where Xq = xo and Xq 

same line of steps as in Section [II. However, as we have seen in Section IV, the levels not satisfying the condition Eq. 
( [l9| ) become higher in energy by O(l) than those which do. This is because the roots of the Bethe ansatz equation, 
xj or xq^j, for j finite, may not remain close to x^ . So, Eq. (CI) is useful only to the levels which do satisfy Eq. 
(ES). Repeating the algebra, the sum Eq. (27) is then found to take the form 



N f dr 
S[f] = 2^if f( x ) R "W— + ( m + W2)/(4) - ( m + ^-/2)/(^o ) 

+/'(x )$ + /"(-to)* + 0(N- 3 ) 
where $ and * are short notations for 

$ = a 2 TTY 2 (y + , y-)/N + a^ 2 Y A {y + , y-)/N 2 , 
*= l-a 2 2 ir 2 Y 4 (y + ,y_)/N 2 



(C2) 

(C3) 
(C4) 



with 



Y2n(y+,y-) = -iJ+(y+ + i - 2mi) + iJ y l l, {y- - i - 2mi) 

fe=i 

- *(! + 2m - 2 Pk)) n - (y- + *(i - 2m - 2 ^-))"} 



(C5) 



fe=i 
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Here, J±(y) are a generalization of Eq. (p9|): 



,<»>,.,, _ 1 r (v + t)"-(y-t)\, ^ i „ !_.,„+, 



When y + = y_ and m = 0, (— l)™l2n(z/, 2/) reduces to Y2„(y) of I. 
Using this general sum formula, we obtain the solution for Rn as 



R N (x) = *.(*) + - + -£ (-T( + - + 0(N~% (C7) 

TV TV *-f 1 + aM V^o / V 2i x a; 2 . / v ; v y 



1~4 

n^O 

and the self-consistency equations for , a2 and 04 as 



± im? + y ± ^- = ±t(l - g)7r + 2H - A + 2 V y^— - 







+(1 + * - y+~y-) In |xo| + (2zto - ^±J^) ^ 
TV 2Ni J 1 ' l y 2 ' TV 

+2e Xo 1 ($-*/so)TV- 1 + 0(TV- 4 ), (C8) 



- a^ 1 - i^fco) + 2reQXQ 1 N~ 1 + 0(TV~ 3 ), 
2a 4 a- 3 = iZ'^(x ) + O(N- 1 ) 



(C9) 
(CIO) 



where Zoo(x) is given in Eq. ( |l7| ) and y+ — ?/_ = 4ir which follows from Eq. (C8) has been used in Eq. (C9). From 
Eq. (|C8|), one gets, after using Eq. (|l§l), 



y± = (±2r + m)i + n^eo^/xo - tf/x§) + 0(TV~ 3 ). 
Applying Eq. @ with Eqs. (jC^) and ( |cTl| ) to Eqs. (§) and (§), we then obtain after some algebra, 
A£ = 2rS'(6) + H"(6)(xo x $ - -V 2 *) - Z"'(b)x^ 2 V + 0(N- 3 ) 

for the chain and 



AF 



2r(V - H(A - &)) - S'(A - fr)^ 1 * 



x,7 2 *) - S"(A - b) X Q 2 * + 0(TV~ 3 ) 



(Gil) 
(C12) 
(C13) 



f or th e lattice model where xo = — exp(A — b) and S(b) is defined in Eq. (|l9|). The 0(1) terms in Eqs. ( C12| ) and 
(|C13|) imply that the critical field V c as a function of 77 is given by S'(6) for the chai n and by 5 (A — b) for the lattice. 
These 0(1) terms are offset by the vertical field at the phase boundary. Eqs. (C12) and ( C13| ) contain, through $, 
an 0(N~ 1 ) term proportional to Y 2 (y± = ±2ri + mi) which is purely imaginary. Thus the real parts of the energy 
gaps at the antiferromagnetic phase boundary are 0(TV~ 2 ). 

Up to here H is general. For simplicity, we now specialize to H = at corresponding v alues of th e crit ical vertical 
field. When H = 0, x Q = -e\ b = and a 2 = e x /e 2 + C^TV^ 1 ). Setting b = in Eqs. (|Cl|) and (|C1^ ), one notes 
that the imaginary 0(N~' 1 ) terms all vanish simplifying the results considerably. Final expression for AE and AF 

then becomes as given in Eqs. ( ff2| ) and (|7C|), respectively, where y[ r ' m ' VS> is the value of Ki(y+,y_) evaluated at the 
zero-th order solution for y±, i.e. Y±(y± = i(±2r + to)). 



APPENDIX D: PARTITION FUNCTIONS AT H = AND V = V c 

In this Appendix, we derive Eq. (|75|). Consider the lattice model on M x TV lattice with M and TV even. Then 
from Eqs. (73) and (|74j), we have 



z F = p xr 



(Dl) 



with yjf'" 1 '^ given in Eq. ( |7l| ) and the sum is for r = 0, 1, . . . under the restriction Eq. (^9|). Considering first (2r + m) 
even sector, from Eq. (71), one sees that a particle (hole) at position pk (h k ) contributes p-( 1 + 2 J') 2 (p( 1 + 2 J') 2 ) with 



1G 



< i < r — m/2 — 1 ( j > r — m/2). Thus the free fermion grand-canonical partition function for the right-hand side 
excitation is 

i — m/2— 1 oo 

ss(*)= n a+^p- (i+2i)2 ) n (i+^ x p (i+2i)2 ) 

3=0 j=r—m/2 
r-m/2— 1 



J=0 



where 



C(z) = + z- 1 p (1+2j)2 ) = ^ C n 2r- n . (D3) 

j=0 n=0 

Thus the coefficient of z° in £3(2:) can be written as C r - m / 2 n^=cT^ 2_1 P _< - 1+2 ^ 2 - The left-hand side excitations give 

the same result with m -> -to. Using the fact that E^o^ 2 "^ 1 + 2j) 2 + YJ'jt™' 2 ' 1 ^ + 2j) 2 = 8r 3 /3 - 2r/3 + 2rm 2 , 

we note that the product in the second line of Eq. ( |D2| ) is cancelled by that coming from the first terms of y^ ,m,V \ 
Thus we have 

pW"'* = C r _ TO/2 C r+m/2 . (D4) 



Summing Eq. (D4) over m from — 2r to +2r in steps of 2 and over r > 0, we find the partition function in the (2r + m) 
even sector as (C*(l) 2 + C(-l) 2 )/2 = Zi(p) + Z 3 (p). 

Next consider the (2r + m) odd sector. In this case the free fermion grand-canonical partition function for the 
right-hand side excitations is 

r-(m+l)/2 00 

2 4 (2)= [] i 1 +^~ if ) II (1 + ^ 1 P 4/ ) (D5) 

j=0 j=r-(m-l)/2 

while that for the left-hand side is 

r+(m-l)/2 00 

S 4 (z)= J] (1 + zp- 4 ^' 2 ) J] (l + ^ X P 4j2 ). (D6) 

3=1 j=r+(m+l)/2 

Note that the j = term is excluded in Eq. ( |D6| ) due to the strict inequality in Eq. ([l9|) . Defining 

00 00 
D(z) = l[(l + z-V j2 ) = J2 D " z ~ n (° 7 ) 

j = 71 = 

00 00 

D(z) = Y[(l + z- 1 p^) = Y / D n z- n , (D8) 



^ P^ 4 = -Dr-(m-l)/2-Dr+(m-l)/2 (D9) 



we can write 

.,.,„ r 

P 

•p 

where E^™ /2 ^ 1/2 ( 2 .?') 2 + Ej=i" /2_1/2 ( 2 i) 2 = 8r3 /3 - 2r/3 + 2rm 2 has been used. We then sum Eq. @ over m 
from — 2r + 1 to 2r + 1 in steps of 2 and finally over r > to obtain the remaining piece of the partition function as 
(D(l)D(l) + D(-l)5(-l))/2 = Z 2 (p). We thus arrive at Eq. ffih for M and TV even. For other parities of M and 
N, we repeat above steps with only minor changes and find Eq. (|75|). 
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FIG. 1. (a) Phase diagram in fields of the asymmetric XXZ chain for A = 2.5. The solid curve is the antiferromagnetic phase 
boundary and the dashed ones are the ferromagnetic phase boundary. The region between the curves is the critical phase. The 
point A (B) is the H = H a , V = (H = 0, V = V c ) phase boundary disussed in Section @. (b) The same as in (a) for the 
isotropic 6- vertex model (F-model). 



FIG. 2. Examples of energy levels for N = 20 and Q = 10 (r = 0). Filled circles denotes integers included in {Ij}. (a) 
The ground state (0,0,0) (b) The 1-shifted level (0,1,0) (c) An excited level from (b) which has two particle-hole pairs at the 
right-end of the Fermi sea (with positions pi=2, P2=3, hi = l and h,2=2) and one particle-hole pair at the left-end (with pi=l 
and hi — 1). 
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Fig. 1(b) 
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